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ABSTRACT 

The recent extreme ultraheavy cosmic-ray observations (Z > TO) arc compared 
with r-process models. A detailed cosmic ruy propagation calculation is used to .runsform 
the calculated source distributions to those observed ut the earth. The r-process produc- 
tion abundances are calculated using different mass formulae and /3-rate formulae; an 
empirical stimate based on the observed solar system abundances is used also. There is 
the continued strong indication of an r-process dominance in the extreme ultra-heavy 
cosmic rays. However it is shown that the observed high actinide/Pt ratio in the cosmic 
rays cannot be fit with the same r-process calculation which also fits the solar system 
material. This result suggests that the cosmic rays probably undergo some preferential 
acceleration in addition to the apparent general enrichment in heavy (r-process) material. 
As estimate also is made of the expected relative abundance of superheavy elements in 
the cosmic rays if the anomolous heavy xenon in carbonaceous chondrites is due to a 
fissioning superheavy element. 
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I. INTRODUCTION 

Experimental measurements of the extreme ultra-heavy cosmic rays (Z > 70) are of 
great interest in the study of nucleosynthesis. Ultra-heavy nuclei are synthesized 
primarily by neutron capture on either a rapid (r) or slow (s) timescale relative to 0 -decay 
rates (Burbidge et al. 1957, Seeger et al . 1965). ^See also the discussion of the generalized 
n-process (Blake and Schramm, 1976) j. By examining the relative abundance of the var- 
ious ultra-heavy nuclei in cosmic rays, one can judge the importance of each process in 
the cosmic-ray source. Conventional wisdom attributes the origin of the bulk of the 
cosmic rays to supernovae and, thus, since r-process nucleosynthesis is believed to occur 
in supernovae whereas s-process nucleosynthesis is not, the ultra-heavy cosmic rays are 
expected to be largely r-process in origin (Schramm 1972, Blake and Schramm 1974, 
Hainebach et al. 1976). The ultra-heavy data currently available (Fowler 1973, Price and 
Shirk 1975, Price 1976, Fowler et al., 1977) definitely show an r-process character in that 
the Pt charge group (74<Z<79) is substantially more abundant than the Pb charge group 
(80 £ Z < 84) and the actinide nuclei are prominent. The publication of the Skylab results 
by Price and Shirk (1975) increased by fifty percent the number of ultra-heavy nuclei 

observed and has the advantage, relative to the other work, of being acquired above the 

2 

atmosphere, albeit inside the Skylab hull ( -1 gm/cm ). 

The nuclear parameters (mass law, -rates, fission rates, fission fragment mass 
distribution) necessary for r-process calculations are not known but must be estimated by 
theoretical analysis and extrapolation from known nuclei. Examination of nearly pure, and 
young r-process material in the ultra-heavy cosmic rays (if this is established 
unquestionably to be the case) can reveal much about the details of the r-process, in par- 
ticular, and supernovae in general. 


1 


I 


l 


I 


I 


I 


In this paper the implications for nucleosynthesis theory are considered of the 
ultra-heavy cosmic-ray experimental results. First the basic uncertainties in the r- 
process calculations are discussed. Next various r-process source-abundance distributions 
are used in a galactic propagation calculation and compared with the Skylab results of 
Price and Shirk (1975). This is the most complete cosmic-ray propagation calculation in 
this mass region to date and includes a complete reaction network. It is found that the 
disparity between the high actinide/platinum ratio in the cosmic ruys and the lower ratio 
in solar system material cannot be explained solely by propagation effects operating on 
the r-process abundances resulting from any single mass law and /?-rate formalism. 
Preferential acceleration effects are required. These are discussed and comment is made 
on uncertainties in the solar system abundances of Pb, Bi and the actinides. The 
importance of improved charge resolution in ultra-heavy cosmic-ray measurements and 
the question of superheavy element abundances are addressed. 
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II. r-PROCESS CALCULATIONS 

Because the nuclei involved in an r-process event are far removed from the valley 
of /J -stability (Burbidge et al. 1957, Seeger et 1965), the relevant nuclear parameters 
are not obtained from experimental measurements. Rather, the mass law and -rates 
used in r-process calculations are estimated by theoretical analysis and extrapolation from 
known nuclei. 

The Myers and Swiatecki (1966) liquid-drop mass law with the modification to the 
shell correction given by Myers and Swiatecki (1967) is used in the r-process calculations 
discussed in this paper. The newer droplet model is available (Myers and Swiatecki 1969, 
Myers 1976) and has been used extensively in calculations also. However the differences 
in the calculated abundances as a result of using the droplet model instead of the liquid 
drop model are small and immaterial for present purposes. The liquid-drop model is used 
in order that a direct comparison can be made with the previous calculations and 
discussions ( viz . Schramm and Fiset 1973, Blake and Schramm 1974). The next neutron 
magic number after N = 126 is taken to be N = 184 and the next proton magic number 
after Z = 82 is taken to be Z_= 114 in the liquid-drop muss law. The values of these next 
magic numbers are not known from experiment and are subject to theoretical debate (cf. 
Nilsson 1972); this is discussed further below. For nuclei in the uctinide region and below, 
the selection of the next proton magic number to be Z = 1 14 or Z = 126 is not significant; 
it is crucial only to a discussion of the r-process synthesis of superheavy elements. 

Two modern descriptions of the /9 -rates are available. One is basically the Fermi 
theory and has been discussed from the point of view of its r-process applications by 
Senbetu (1973). The second /9-rate formalism is the Gross theory of /9 -decay (Takahushi 
and Yamada 1969, Kodama et al. 1970, Takahashi, 1971) and has been described and used 
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in r-process calculations by Kodama and Takahashi (1975). The Senbetu formalism is a 
"simple" description which requires an arbitrary log ft value and averages overall stutes 
without special weighting. The Gross theory includes more of the relevant physics. 
Because nuclei far from 0 stability have large Q values, their decays proceed to many 
final states. The Gross theory was devised to treat this situation and, instead of dealing 
with individual transitions, it directly treats the 0 strength functions. However, it must 
be emphasized that the nuclei of interest in an r-process are far removed from the region 
where experimental data presently exists. Any 0 -rate formalism thus is a large extra- 
polation into the unknown. It turns out that the major difference in the re 1 ults of the 
present calculations between the two different 0 -rates (for r-process abundances of 
interest 'n the present paper) is that the abundance of the actinides relative to the plati- 
num peak is a factor of — J larger using the Gross theory. 

In the Senbetu /3-rate formalism, the transition probabilities, leg ft, are chosen to 
be some average value for all transitions. The value chosen (usually between 5 and 6.5) 
does not in any way affect equilibrium r-process calculations since it treats all transitions 
equally. It is of significance in dynamic calculations (Schramm 1973, Hillebrandt and 
Takahashi 1976) but even then has little effect on relative abundances in a given mass 
region. It is possible to have the log ft value vary with A when using the Senbetu 
formalism (Blake ar.d Schramm 1974). An examination of the 0 -rates of nuclei near the 
valley of 0 stability indicates that Aar80 may be a transition point from log ft ar 5 to log 
ft at 6.5. However, far from the valley of stability, there is no experimental data 
available to guide the selection. In the case of the Gross theory of beta decay, the log ft 
values are not free parameters. 

As indicated above, variations in the mass law for a given description of the 0 - 
rctes also strongly affects the r-process abundance predictions (Seeger 1967, Mathews and 
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Violn 1976). A marked effect '(or present purposes) arises If the next neutron magic 
number after N = 126 Is N = 164 rather than N * 184; such is the case in the mass law of 
Seeger and Howard (1975). A magic 'umber at N = 164 yields an abundance peak on the 
valley of & -stability at A a 250, in the region of the heavy actinides. The result is that 
the heavy actinides would be produced in much greater abundance than the low -mass 
actinides, the short-lived trans-bismuth nuclei, lead, and bismuth. Mathews and Viola 
(1976) show the shape of an r-process production curve using the Seeger and Howard (1975, 
1976) mass law. Such a production curve is not a good fit to the empirical data. 
Furthermore, theoretical opinion appears to strongly favor M = 184 as the r-»xt magic 
number after N = 126. Mathews and Viola (1976) also give an r-process production curve 
using the empirical mass law of Viola et al. (1974); this mass law gives a production in the 
actinide region within a factor of two of the platinum peak. The empirical evidence 
argues strongly aguinst this mass luw; actinides would have a much greater abundance 
relative to the platinum peak than is observed. 

The results presented in this paper are from static (constant temperature, constant 
neutron density) r-process calculations. Because the astrophysical setting of the synthesis 
of the neutron-rich nuclei attributed to r-process is unknown (cf. Schramm and Norman 
1976), a more complex (dynamic) calculation simply adds additional free parameters. 
Although the dynamic calculations are more "realistic” and might eventually permit 
identification of the r-process site, such calculations do not give a convenient focus on the 
undei lying nuclear physics because of the substantial additional complication caused by 
the hydrodynamic parameterization. It has been shown (Schramm 1973) that the general 
character of the smoothed dynamic solution is similar to that obtained with a 
straightforward static calculation in which the temperature and neutron density are 
chosen in such a way as to produce the observed abundance peak locations. A good fit 
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orcured in the present calculations for * 1.8 arul log J() neutron density * 28. In a 
dynamic calculation the neutron flux decreases with time (Schramm 1973, llillebrandt and 
Takahashi 1976, Schramm and Norman 1976). However, if all the actinides are 
synthesized in a given event, whether the r-process path terminates due to neutron in- 
duced fission or neutron exhaustion is irrelevant to the relative actinide abundances. It is 
possible for un r-proeers event to tern. mate due to neutron exhaustion iust as »he 
actinides begin to be produced and thus yield, for example, a Th/Pt ratio of arbitrary 
value between the standard cyclic solution value and zero. However, the observational 
data gives such ratios as Pu/Th, U/Th, actinides/Pt, aetinides/Bi which suggest that the 
relative r-process production abundances in the actinide region are quite uniform. This 
would not be the case if the r-process ♦erminates just as the actinide regior is reached. 

In addition to th- v ove theoretical calculations, a set of initial r-process 

abundances implied the present solar system r-process abundances have been 

•>()•( 

generated. Specifically, these are derived from the present-day abundances of “ Bi 

(from the decay of “ 37 Np), ^“Th, ‘' 33 U, and “ 38 U using a constant rate-of-synthesis 

model for Galactic chemical evolution, and from an assumption of equal initial abundances 

232 238 

of all even r iclei (the Th and U chains), and equal initial abundances of ull odd 
nuclei (the ^ 3 Np and , ‘ >35 U chains). 


J 
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III. PROPAGATION CALCULATIONS 


It is necessary to account for alteration of the competition of the cosmic rays due 

to collisions with ambient particles in the interstellar medium and due to r adioactive 

decay. The present calculations depend largely on the semi-empiricul formulae for 

nucleur breakup react ; on cross sections given by Silberberg and Tsao ( 1973a, b). However, 

for fission and spalKtior. of trans-bismuth nuclei, the fission-corrected semi-empirical 

formulae of Schramm (1972) have been used. For peripheral reactions on trans-bismuth 
238 

nuclei other than U, the cross sections used are those of the corresponding reactions 
238 

for U, but scaled by the relevant total inelastic cross sections and fissilities. Total 
inelastic loss cross sections were taken from Kirkby and Link G966). 

A sieady-state model of cosmic ray transport through exponentially-weighted slub 
distribuitons wus used for the propagation calculations discussed in this paper. Only 
energy independent exponential path length distributions were considered and no attempt 
wus made to reproduce the secondary to primary ral.'*. declining with increasing energy, 
which have been observed in the light, medium, and heavy cosmic rays (Juliusson 1974). 
The effects of energy jss are taken into account by assuming stopping powers appropriate 
to a neutral gas with solar system compositon. The resulting model is equivalent to the 
leaky box model of Gloeckler and Jokipii (1969) with a leakage time for particle: from the 
Galaxy inversely proportional to velocity. 

A total of 177 species with A > 170 were included in the present calculations. 
Nuclei with alpha, beta, and fission decay modes were included if T J/2 > 100 yr. Pure 
electron-capture nuclei were included if T^fEC) > 10 h and T^l /?*) £ 100 yi. The 
cross sections needed to treat electron capture nuclei were calculated using the Brinkman 
and Kramers (1930) formula for non-radiative capture and the formula given by Bohr 
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(1948) for electron stripping. Or.e-electron capture cross-sections in the charge region of 
interest in this paper are large at all energies likely to be of interest in the near future. 
Thus most of the electron capture species in this mass region should decay during 
propagation unless they have half-lives comparable to or longer thun the age of the cosmic 
rnys. This happenstance greatly reduces the number of electron capture nuclei that need 
to oe included in a detailed study of cosmic-ray propagation in the Pt region and above. 

The cross sections of all possible reaction channels were added together in es- 
timating the interaction length for one species going tc unother species. The mean inter- 
stellar density enters the calculation in relating decay half-lives to interaction lengths in 
the interstellar gas. One of the unique features of the present calculations! method lies in 
the diagonalization of the interaction matrix. This allows an accurate treatment of fast 
decays and all higher order production modes within the puthlength step chosen for the 
integration of the partial differential equations describing production and loss of each 
nucleus along its corresponding energy-loss characteristic. This diagonalization procedure 
is relatively simple to carry out because the interaction matrix is almost upper triangular 
when the nuclei are ordered in terms of increasing A and decreasing Z. This allows the 

eigcnrow and column vectors .o be calculated by straightforward Gaussian elimination. 

« 6 

The source spectra used in this calculation were of the form Q.(E ♦ E (( ) where 
Q ( is the relative abundance of species i; E is kinetic energy/nucleon; and E fl = 400 
MeV/nucleon. This form of the source spectra, when modulated, produces a good fit to 
the observed energy spectra of light and medium cosmic ray nuclei (Garcia-Munoz et al. 
1975b). 

The recent studies of the age of the cosmic rays, using 10 Be (Garcia-Munoz et al. 
1975a, 1977a) led to the parameters used in the present propagation calculations: an 

exponential path-length distribution with a 5 g/cm (of hydrogen) leakage mean-free-path 
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(6 g/cm including interstellar helium). Cosmic ray propagation calculations frequenMy 

.3 

assume an interstellar matter density of at l atom/cm . However, Garciu-Munoz et al. 

.3 

(1975a, 1977a) found a much lower value, ar0.2 atoms/cm , giving a leaknge lifetime of 2 

7 -3 

x 10 yr. We have done propagation calculations using both 1 atom/cm and 0.1 

3 

atom/cm and present the latter results in this paper as being more consistent with the 
results of Garcia-Mur.oz et al. (1975a, 1977a). This lower density leads to a cosmic ray 
leakage age of 3.2 x 10 yr for a velocity of /9* 1. 
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IV. RESULTS 


•***>*0 


*AC£ 


Not 

The available experimental data and the results of the propagatior - cafKulation done 




with several initial relative ubundances are shown in Table 1. The statement that the 
experimental data strongly favors an r-process source can be seen explicitly. The high 
actinide relative abundance and the presence of trans-uranics cannot result from the 
acceleration of solar-system material. Furthermore, it can be seen that the Gross theory 
/? -rates yield a reasonable fit to the measured U/Pt and Trans-U/Pt ratios, especially to 
the Skylab data. The Senbetu yj-rates yield U/Pt and Trans-U/Pt ratios approximately a 
factor of three lower than the Skylab results. 


oio 91* 0*10 

In Table 2, the solar system abundances of Th, U, and U (Cameron 1973) 

are presented along with the calculated r-process abundances of these isotopes using the 

two /J-rate formalisms. In order to test the consistency of the calculated abundances of 
232 235 238 

Th, U, and U, an extreme model for galactic evolution was taken in which it was 

a 

assumed that all r-process elements were created in one burst of nucleosynthesis 20 x 10 y 

ago, i.e. the Hubble age for the Universe. This assumption assures that the maximum 

amount of actinides will have decayed, and gives the lowest possible actinides/Pt ratios at 

uie time of formation of the solar system. One sees that the aged Senbetu abundances are 

k smaller than Cameron's values for solar system abundances, indicating that it is necessary 

to construct a less extreme galactic evolution model which agrees with tne presently 

observed abundances of these elements (cf. Hainebach and Schramm 1977). However, 

232 

after aging the Gross-theory abundances, the resulting Th abundance is still a factor of 

~2 higher than its solar system abundance. In addition, the Gross theory abundance of the 

209 209 

Bi decay chain is also a factor of — 3 higher than the solar system abundance of “ Bi. 


-11- 


I 




u 

-2 

< 

H 



‘O' 

o> « ^ a> 

o o o o o 

• • • • • 

o o o o o 


m 

«* 

o o 

o o 

O O <-5 'J- 

o o o o 

• • • • 

o o o o 







co *r 



CO 

co 


• 

o 


© © 

c- N 
• • 

cm uo 
• • 

© 

• 

® 

£ 

44 



o o 







*-« CO 

, W' 

W W 

wmm 

CM 


*r 

<e 

CM CM 

uo uo 

00 o 

CM 

CM 

a. 

o 

o 

• • 
o o 

• • 
o o 

• • 
o ^ 

• 

o 

© 




-12- 


t t 


-Process Actinide Abundances 


<2 


CO 

in 


E 

u. 

t- 

%n 

cm 


s 

* 

V 

<0. 

2? 

8 


o> 


6 SfS 


in 

o 

<r> 


ao 

T 

o ^ 

s - 2 

x © 

° CM © 


o 

© 

«*■ 

CM 


E 

O 

1/3 

•r 

CM 


CM 


I 


g . £■ 

r 8 1 

- i 5 

•2 c ft 

•- 2 © 

— co 6b 


in 

© 

V© 

© 


CO 03 

CM ** 


o. 

z 


co 

CM 


3 

CO 

CO 

CM 


E 

u* 

UD 

in 

CM 


8 

8; 





xz 

F 

Cf, 

Fm 

w 

03 





03 

CM 

in 

S3 

o 





CM 

in 

in 

V 






CM 

CM 

CM 

Ta 

X 



00 





3 

in 



1 

CM 

CO 


• 


0) 


o 

in 

o 

S“ 

E 

<•** 

■ C 



CM 


o 

c 

o 

o 

c 

$ 


o 

X 

© 

in 

00 

•■w 

3 

bfl 

03 

• 

o 

• 

o 

00 

d 

© 

CM 

CM 

*r 

CM 

in 

CM 


8 


v x 

l «V 

«?s 

~ 3 C 

'.SScB 


CO 

in 


© 

o 


co 

co 

o 


u- 


CM 

CM 


< 

C- 

CM 


3 

C- 


3 

x 

o 

CM 


E 

U 


E 

«£ 

CO 

cm 


E 

O 


cc 

*T 

CM 


J- </J 
OC 0) 

s & i 

e « , g 

E 3 % 

3 *3 A 
cn •- « 


CO 


00 

in 

o 


CO 

<£> 


CM 


o 

o- 

co 

CM 


X 

03 

o 

CM 


<r 

co 

CM 


CM 

CO 

CM 


3 

X 

03 

CO 

CM 


X 

in 

co 

CM 


3 

X 

CM 

CM 


X 

00 

CO 

CM 


8 c 
*1 

X 

03 

CM 

D 

in 

X 

00 

Ln x: 

o 

CO 

co 

CO 

X o 

CM 

CM 

CM 

CM 


3 

x: 

o 

0) 

£ 


X 

03 

o 

M 


£Z 

• • 

• 

F 

D 


CM 

m 

00 

CO 

CO 

CO 

CM 

CM 

CM 


- 13 - 


! 


Some chains may go higher, cf. Schramm and Fowler (1971). 

t Cameron (1973): Si = 10 6 ; values at the time the solar system condensed 4.6 x 10 y yr ago. 

* Normalized to 195 Pt - 0.473. the solar-system abundance. 

*15 x 10® yrs before the solar system condensed, i.e.. ~20 x Id® yr before the present, «hieh is the Hubble K* 



Thus, it is not possibl* to construct a model of galactic nucleosynthesis using the dross 

theory of 0 -decay which is consistent with the Cameron (1973) solar system abundances 

°09 232 238 

without assuming that both “ Bi and Th are depleted relative to U by chemical 

fractionation during the formation of tne solar system. 

An examination of the meteoritic data collected in Mason (1971) reveals that it is 

**32 

unlikely although not inconceivable that the Th abundance adopted by Cameron (1973) 
is in error by as much as a factor of 2. The solar “Th abundance was recently reduced 
by a factor of four and r.ow agrees with the present meteoritic abundance as a result of a 
new value for the Th oscillator strength derived from laboratory measurements (Anderson 
and Petkov 1975). The meteoritic abundances of Bi and U are better determined. Both 
elements have abundance variations of less than 10% in Cl chondrites (cf., Krflhenbuhl et 
a]. 1973). Bi is more volatile than Th or U as evidenced by the fact that it is depleted in 
C2 chondrites by a factor of about 2. However, the degree of chemical fractionation of Bi 
and Th required to make the level of actinides produced in the Gross-theory calculation 
consistent with the solar system abundances of these elements is unlikely. 

In order to explain the U/Pt and Trans-U/Pt ratios observed in the ultra-heavy 
cosmic rays, one needs a total actinide abundance similar to that predicted by the present 
static r-process calculation using the Gross-theory 0 -rates (see Figure 1). Increasing the 

odd to even A abundances in the actinide region while maintaining the same total actinide 

232 209 

abundance eases the Th problem but makes the required fractionation of Bi even 

232 238 

greater. It is possible to reduce the Th/ U production ratio about a factor 2 to 

209 

obtain agreement with the solar system ratio but one is still left with too much Bi. 
The basic problem is that there are more long-lived actinides with odd A then even A 
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Fig. 1. The calculated r-process abundances for A >180 are shown using both the 
Senbetu /? -rates and the Gross theory /?-rates. The liquid drop mass model was used and 
a temperature of Tg = 1.8 and log neutron density = 28 (which fits the observed solar- 
system peak platinum at A = 195). Adjacent values of A are connected by lines to guide 
the eye except at high atomic weight where spontaneous fission (on the valley of stability) 
eliminates some of the adjacent mass numbers. 
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contributing to the abundances of the ultra-heavy cosmic rays, and for published estimates 
of odd to even production ratios it is not chronologically possible to increase the 
production abundance of actinides sufficiently to explain the cosmic ray abundances 
without o er-prcducing Bi (further details on the chronological restrictions on r- 
process abundances are to be discussed in Anglin 1977). Thus, unless the cosmic ray data 
is systematically in error, it is not possible to explain the high actinide abundance in the 
cosmic rays and solar system actinide abundances with the same r-process. One must look 
for another explanation of the high actinide abundance of cosmic rays. 

One possibility is that selection effects, which depend upon the atomic properties 
of the elements, operate at the source prior to and/or during acceleration. An argument 
against such selection effects operating on ultra-heavy cosmic rays is that while they 
appear to occur for elements lighter than Si, the cosmic ray source abundances of Si, Ca, 
Fe, and Ni are remarkably similar to those found in the solar system (Garcia-Munoz et al . 
1977b). However, selection effects can't be ruled out as a possibility at the present time. 
Kristiansson (1971, 1972, 1974), Casse and Goret (1973), Havnes (1973), and Casse et al . 
(1975) have considered various selection effects based upon nuclear charge Z, electron- 
impact ionization cross-section a , and the first ionization potential I. The most recent 
paper, Casse et al . (1975), examines all three possibilities. Their analysis suggests that 
element selection dependent on ionization potential is the most reasonable, with 
overabundance increasing with decreasing ionization potential. Values for I for the 
elements of major interest in the ultra-heavy cosmic rays data discussed here are given in 
Table 3. 

It is clear from the run of values in Table 3 that the actinides would be favored 
relative to the elements in the platinum peak if the abundances huve the proposed de- 
pendence upon ionization potential. In order to quantify the result somewhat, the data 
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TABLE 3 

Ionization Potential of Selected Elements 


Element 

Ionization Pol 

Os 

8.70 

Pt 

9.00 

Au 

9.22 

Pb 

7.42 

Th 

7.50 

U 

6.19 

Pu 

5.71 


From the compilation of Kaye and Laby (1973) . 
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presented by Casse et al. (1975), (their Figure 1-e relating relative abundance to ioni; . 
tion potential) have been fit with a straight line (on a semi-log plot) for 5 eV£ I <15 eV. 
(The data for H, Ne, and He were not used.) The straight-line fit predicts the following 
enhancements for U and Pb relative to Pt: 

U/Pt enhanced by a factor of 2.5, 

Pb/Pt enhanced by a factor 1.5. 

It is interesting to note that if such ionization-potential-dependent enhancements 
actually occur, using the results of either the r-process calculations with the Senbetu 0 - 
rates or the implied solar system r-process abundances after propagation yield a good fit 
to the cosmic ray (except perhaps Pb/Pt, but see below). However, the Gross theory 0 - 
rates now yield too high an actinide/Pt ratio. In this regard, note that the solar system 
abundances of stable r-process elements are fit better using the Senbetu 0 -rates than the 
Gross theory /J-rates (Blake and Schramm, 1974). It is clear ti nt an understanding of 
preferential acceleration or other selection effects is crucial to further progress in 
relating observations to theoretical r-process calculations. 

Note that the Pb/Pt ratio would be increased by preferential acceleration as dis- 
cussed above. Preferential acceleration thus would not be able to alter a solar system 
composition source into an enriched Pt peak relative to Pb, but in fact would deplete Pt 
relative to Pb. One further point to remember is that preferential acceleration of the 
radioactive nuclei between Bi and Th would further increase the Pb/Pt ratio since these 
radioactive nuclei will decay to the Pb peak, thus bringing the Senbetu # -rates r-process 
abundances, with preferential acceleration, into better agreement with the Skylab data. 
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The possibility that r-process events differ in their production ratios and that the 
cosmic ray source and solar system material each comprise different distributions of r- 
process events should be considered. However, the sharpness of the solar system r-process 
peaks, e.g. the Pt peak, argues against this: if the location in mass number of the 

peaks from various events differed much, the composite peak would be much broader than 
observed. One might instead propose tha' the cosmic ray source consists of a few peculiar 
supernovae ( e.g. Chevalier's suggestion, discussed below) which contribute little to 
Galactic chemical evolution, but in which the r-process environment persists longer and 
puts more material into heavier nuclides, like the actinides, than does the more common 
(in this model) r-process. This rare r-process would be the ordinary cyclic solution r- 
process (Seeger et al. 1965) in which material eventually reaches some value of A where 
neutron-induced fission occurs, the fission fragments theieafler serving as r-process seed 
nuclei. In order to account for the high cosmic ray actinide/Pt ratio, it would have to be 
assumed that the usual theoretical cyclic r-process is actually rare in nature and that 
solar systen. material came from a non-cyclic r-process, which gave a lower actinide/Pt 
ratio. But there is evidence that solar system material is the result of a cyclic r-process, 
viz , the abundance "hump" in the rare earth region which could be the result of fission of 
material at the N=184 magic number (Schramm and Fowler 1971). Thus while a non- 
standard r-process might explain (in an ad hoc fashion) the high cosmic ray actinide/Pt 
ratio, such an explanation, upon examination, is not very promising. 

It is possible also that uncertainties in the measurements might have significantly 
affected the published results. Price and Shirk (1975) give a median charge resolution of 
AZ = 2.5. They take Au (Z = 79) as the upper bound of the platinum peak and Hg (A = 80) 
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is the lower bound of the lead peak. The r-proeess is the major source of Au and even of 
Hg. Thus substantial amounts of r-pri.-’ess material may well be assigned to the lead peak; 
i 9 also Pt itself has Z * 78. The lead peak is overwhelmingly Pb (Z = 82) and Bl (Z = 83); 
thus little actual lead-peak (Le., s-process) material would be expected to be misussigned 
to the platinum peak. 

The r-process indirectly produces substantial amounts of leud-peak material. As 
;ne actinide production increases, the lead peuk does also, due to the fact that the lead 
peak (as produced In the r-process) largely arises from the rapid decay of . uclei with 
84SZS89, and not from the direct production of Pb and Bl. The dilemma, that a 
substantial lead peak can be due either to enhanced actinide abundances or to the 
presence of s-process nuclei, can be experimentally addressed with a good charge- 
resolution measurement (UK-6 and HEAO-C) able to separate Pb from Bi. For nuclei 

7 

created in the r-process <- 10 yr prior to observation, Bi/Pb^0.5, whereas in the solar 
system (Cameron 1973), Bi/Pb*0.035. 

Some s-process material, originally present in the envelope of supernovae, would be 
expected in the ultra-heavy cosmic rays, but the ratio r/s is expected to be of the order of 
10 or 20 (Hainebach et ah 1976). Chevalier (1976) has suggested that heavy element 
acceleration occurs in faint supernovae which result when a star, wh. v n would have been a 
normal type II supernovae, loces its envelope before exploding. If such is the case, the r/s 
ratio might be even larger. 

o 

As mentioned in Section III, the leakage mean free path of 5g/cm of interstellar 

3 

hydrogen and interst“'lar density of 0.1 atoms/cm chosen for our calculations lead to a 

n 

cosmic ray leakage lifetime of 3.2 x 10 yr for light nuclei. The actinides however (as 

10 2 
opposed to light nuclei like Be) have a destruction mean free path of ~ 1 gm/cm , and 
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therefore a smaller mean age, perhaps 5 x 10 yr. For this reason, the actinides are not 

particularly useful cosmic-ray chronometers in that they do not give the leakage lifetime 

3 

of the bulk of the cosmic rays. If the interstellar density were ns high as 1 atom/cm , the 
actinide lifetime would be very short indeed. 

There is evidence that a truncated exponential path length distribution provides a 
better fit to the abundances of the secondaries from l.i to Mn than does a pure exponential 
(Garcia-Munoz et nl ., 1977b). The effect of this change is to increase the mean age of the 
observed cosmic rays and the relative abundance of secondaries to primaries given a fixed 
average puth-length. In the ultra-heavy cosmic rays, such a path-length distribution would 
result in a lower actinide/Pt ratio than otherwise end make the observed high ratio even 
more outstanding in the source. 
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V. SUPF.RHEAVY ELEMENTS 

The ultra-heavy cosmic ruy experiments have detected no superheuvy elements 
(Fowler 1973; Price and Shirk 19 V o). However, these limits on the superheavy flux are not 
stringent. 

Recently Anders and colleagues (Lewis et al. 1975; Anders et al. 1975) have 
carried out a beautiful set of experiments studying the rare gases, in particular the excess 
heavy x*non isotopes, in separate fractions of the Allende ..ieteorite. They argue that the 
excess xenon may be the product of a fissioning heavy nucleus and, since this putative pro- 
genitor is shown to be a volatile, actinides are excluded. The results of the detailed 
analysis of Anders et al. (1975) suggest that the superheavy progenitor had a nuclear 
charge of 113, 114 or 115 and that the xenon cannot be due to the presence of primordial 
grains. If the xenon is indeed aue to a fissioning nucleus, the halflife must be such that a 
significant amount could survive the time interval between the last addition of newly 
synthesized nuclei to the solar system material, and the onset of xenon retention. Cos- 

Q 

mochronological studies show this interval, A , to be ~ 2 x 10 yr (Schramm 1974). 

The observations of Lewis et al. (1975) and Anders et al. (1975), and the assumption 

that the xenon is due to a superheavy progenitor, may be used to estimate the expected 

flux of superheavy cosmic rays relative to the actlr ide cosmic rays. At the time of xenon 

238 

retention, the abundance of the superheavy (Sh) relative to U was (Anders et al . 1975) 

r 3. ~ 6 x 10' 4 . (1) 

” *®U 

It is also known fr^m other meteoritic studies (cf. Lewis et a». 1975 and earlier 
refs, therein) that 

fKUJtUlNU PAGE BEANS NOT FIEMEU 
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( 2 ) 


244 

238 


Pu 

U 


*3 x 10 


g 

at the time of xenon retention. Since the value of A is ~ 2 x 10 yr., one has, from 
equations (1) and (2), and the usual prescription for radioactive decay: 


1.8 x 10 * = 


-2 Sh 


7TJ 


Pu I Xe retention 


(3) 


Sh 


7TT 


Pu (production 


exp -2 x 10 8 (I 


1 


Sh 1.2x10 


r 


8 244 

where 1.2 x 10 is the mean life of Pu, and 

superheavy progenitor of the xenon. 


r Sh is the mean life of the fissioning 


Studies of the nucleosynthesis of superheavies in the r-process indicates that the 
ratio {production * s cerlu ‘ n *y not g rea ter than unity (Schramm and Fiset 1973). Using 
this upper limit for the production ratio, equation (3) yields (cf. Schramm 1971) 


Tg h >3.5 x IQ 7 yr. (4) 

If the ultra-heavy cosmic ray age is as large as the 10 Be age (Garcia-Munoz e t ai . 1975a), 

n 

viz . 2 x 10 yr, the calculated r-process abundances of Blake and Schramm (1974) plus 
expression (4) give 


Sh : Trans-U : U = ( <0.06) : 0.45 : 1. 


(5) 
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Many unsuccessful attempts have been made to find superheuvy elements in solar 

system material (cf. Hermann 1974). The negative result is taken to indicate that the 

235 

mean life of the superheavy progenitor of the fission gas is not greater than that of U, 

Q 

i.e., 1.02 x 10 yr. Using this upper limit yield, from equation (3): 


> 4.1 x to -3 . (6) 

production 

, 7 

Again, taking the ultra-heavy cosmic-ray age as 2 x 10 yr gives 

Sh : Trans-U : U = (>0.0004) : 0.45 : 1 (7) 

However, the age of the ultraheuvy cosmic rays may be much less than the l0 Be cosmic 
ray age because of the large cross section for fission during propagation. If the ultraheavy 
cosmic ray age were only 5 x i0 yr, as suggested in Section IV, the abundance ratios, and 
upper and lower limits on superheavies become: 

Sh : Trans-U : U = (0.0005 - 0.10) : 0.75 : 1. (8' 

0 

If the ultraheavy cosmic ray age were as little as 10 yr, this becomes: 

Sh : Trans-U : U = (0.0005 - 0.12) : 1.2 : 1. (9) 


These results may be summarized as follows: Using the Allende analysis of Anders 
et al. (1975) plus reasonable but broad limits on the possible production abundance and 
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lifetime of superheavy nuclei (Schramm and Fowler 1971, Schramm and Fiset 1973) Rives a 
predicted superheavy element flux in the cosmic rays o 


0.0004 S <0.12. (10) 

Recently Flerov (1977) has summarized the results of searches for superheavy 

elements in carbonaceous chondrites and water samples from hot springs. Multiple 

neutron emission has been detected which cannot be explained by the spontaneous fissio-i 

of the U known to be in the sample ( U is the only naturally occurring nucleus with a 

significant probability of spontaneous fission decay). These results may be evidence for 

T 235 1 

the presence of a long-lived r j^3h) ~ T l/2 ( ‘ U) superheavy nucleus with a concentra- 

-14 

tion in the meteorite of the order of 10 g/g. Of course these results need to be 
verified. However one can say already that if superheavy elements with lifetimes near 
that of U are found to exist, and these nuclei are nucleosynthetically produced in 
conjunction with the actinides, then superheavy nuclei eventually should be found in the 
cosmic rays. 
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VI. CONCLUSIONS 

An analysis of the ultraheavy eosmic-ray data, especially those collected during the 
Skylab mission, yields the following results: 

1. For the nuclei in the platinum peak and above, a basically r-proeess source is 
required as noted earlier (Fowler 1973, Price and Shirk 1975). 

2. The disparity between the high actinide/platinum ratio in the (propagated) 
cosmic rays and the lower ratio observed in solar system material cannot be explained by 
the use of a single /9-rate formalism and mass law in an r-process calculation unless 
another process, such as preferential acceleration, enhances this ratio in the cosmic 
ray source. 

3. Using reasonable but broad limits on the possible production abundance and 
lifetime of superheavy nuclei and results of the analysis of the Allende meteorite by 
Anders et al. (1975) gives a predicted superheavy element flux in the cosmic rays of 
0.0004 < ^ < 0.12. 
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